Nonlinear Dynamics of A Damped Magnetic Oscillator 
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We consider a damped magnetic oscillator, consisting of a permanent magnet in a periodically os- 
cillating magnetic field. A detailed investigation of the dynamics of this dissipative magnetic system 
is made by varying the field amplitude A. As A is increased, the damped magnetic oscillator, albeit 
simple looking, exhibits rich dynamical behaviors such as symmetry-breaking pitchfork bifurcations, 
period-doubling transitions to chaos, symmetry-restoring attractor-merging crises, and saddle-node 
bifurcations giving rise to new periodic attractors. Besides these familiar behaviors, a cascade of 
"resurrections" (i.e., an infinite sequence of alternating restabilizations and destabilizations) of the 
stationary points also occurs. It is found that the stationary points restabilize (destabilize) through 
alternating subcritical (supercritical) period-doubling and pitchfork bifurcations. We also discuss 
the critical behaviors in the period-doubling cascades. 

PACS numbers; 05.45.-a, 05.45.Ac 



I. INTRODUCTION 



We consider a permanent magnet of dipole moment 
m placed in a periodically oscillating magnetic field. 
This magnetic oscillator (MO) can be described by a 
second-order non-autonomous ordinary differential equa- 
tion §§, 



16 + b9 + m.Bo cos cji sin 6* == 0, 



(1) 



where the overdot denotes the differentiation with respect 
to time, 9 is the angle between the permanent magnet 
and the magnetic field, / is the moment of inertia about 
a rotation axis, b is the damping parameter, and Bq and 
to are the amplitude and frequency of the periodically 
oscillating magnetic field, respectively. 

Making the normalization ujt -^ 2Tr{t + i) and 6 -^ 
2ttx, we obtain a dimensionless form of Eq. (m), 



Tx — A cos 27rt sin 2'kx = 0, 



(2) 



where 2: is a normalized angle with mod 1, F = 2TTb/lLu 
and A — 27rmBQ/Iuj^. Note also that Eq. (0) describes 
the motion of a particle in a standing wave field [Q-£|. 
For the conservative case of F = 0, the Hamiltonian 
system exhibits period-doubling bifurcations and large- 
scale stochasticity as the normalized field amplitude A 
is increased, which have been found both experimentally 
m^ ^-iid theoretically [0-p[. Here we are interested in 
the damped case of F / and make a detailed investi- 
gation of the dynamical behaviors of the damped MO by 
varying the amplitude A. 

This paper is organized as follows. We first discuss 
bifurcations associated with stability of periodic orbits 



and Lyapunov exponents in the damped MO in Sec. O. 
With increasing A up to sufhcently large values, dynam- 
ical behaviors of the damped MO are then investigated 
in Sec. 111. This very simple-looking damped MO shows 



a richness in its dynamical behaviors. As A is increased, 
breakdown of symmetries via pitchfork bifurcations [^, 
period-doubling transitions to chaos 01, restoration of 
symmetries via attractor-merging crises |I0| , birth of new 
periodic attractors through saddle-node bifurcations [|j, 
and so on are numerically found. In addition to these 
familiar behaviors, the stationary points exhibit a cas- 
cade of "resurrections" uM (i.e., they restabilize after 
their instability, destabilize again, and so forth ad in- 
finitum). It is found that the restabilizations (desta- 
bilizations) occur via alternating subcritical (supercrit- 
ical) period-doubling and pitchfork bifurcations. An infi- 
nite sequence of period-doubling bifurcations, leading to 
chaos, also follows each destabilization of the stationary 
points. In Sec. IV, we also study the critical scaling be- 



haviors in the period-doubling cascades. It is found that 
the critical behaviors are the same as those for the one- 
dimensional (ID) maps ||9|]. Finally, a summary is given 
in Sec. 



II. STABILITY, BIFURCATIONS AND 
LYAPUNOV EXPONENTS 



In this section we first discuss stability of periodic or- 
bits in the damped MO, using the Floquet theory [ p^ . 
Bifurcations associated with the stability and Lyapunov 
exponents are then discussed. 

The second-order ordinary differential equation (0) 
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is reduced to two first-order ordinary differential equa- 
tions: 



y — — Fy + A cos 27rt sin 27ra;. 



(3a) 
(3b) 



Tliese equations have two symmetries Si and 5*2, be- 
cause the transformations 



Si : X 
S2 : X 



X ± -, y 
2 '^ 

-X, y -^ 



-^ y, t 
-y, t- 



2' 



i, 



(4) 
(5) 



leave Eq. (||) invariant. The transformations in Eqs. (Q) 
and (ph are just the shift in both x and t and the (space) 
inversion, respectively. Hereafter, we will call 5*1 and S2 
the shift and inversion symmetries, respectively. If an or- 
bit z{t)[= {x{y),y{t))] is invariant under Si {i — 1,2), it 
is called an Sj-symmetric orbit. Otherwise, it is called an 
S'j-asymmetric orbit and has its "conjugate" orbit Siz{t). 

The phase space of the damped MO is three dimen- 
sional with the coordinates x, y, and t. Since the damped 
MO is periodic in i, it is convenient to regard time as a 
circular coordinate (with mod 1) in the phase space. We 
then consider the surface of section, the x — y plane at 
interger times (i.e., t = m, m: integer). The phase-space 
trajectiory intersects this plane in a sequence of points. 
This sequence of points corresponds to a mapping on the 
plane. This map plot of an initial point zq [— {xq, j/o)] can 
be conveniently generated by sampling the orbit points 
Zm at the discrete time t = m. We call the transformation 
Zm -^ Zm+i the Poincare map and write Zm+i — P{zm)- 

The linear stability of a g-periodic orbit of P such that 
P*(zo) = Zq is determined from the linearized-map ma- 
trix DP'^{zo) of P"* at an orbit point Zg. Here P'' means 
the q-times iterated map. Using the Floquet theory, the 
matrix M (= DP'') can be obtained by integrating the 
linearized equations for small displacements. 



5x 
Sy 



Sy, 
-rsy- 



2ttA cos 2TTt COS 27ra; Sx 



(6a) 
(6b) 



with two initial displacements ((5a;, Sy) = (1, 0) and (0, 1) 
over the period q. The eigenvalues, Ai and A2, of M are 
called the Floquet (stability) multipliers, characterizing 
the orbit stability. By using the Liouville's formula 0, 
we obtain the determinant of M (det M), 



detM 



^-r? 



(7) 



Hence the pair of Floquet multipliers of a periodic orbit 
lies either on the circle of radius e"'"*/^ or on the real axis 
in the complex plane. The periodic orbit is stable when 
both multipliers lie inside the unit circle. We first note 
that they never cross the unit circle, except at the real 
axis and hence Hopf bifurcations do not occur. Conse- 
quently, it can lose its stability only when a Floquet mul- 
tiplier decreases (increases) through —1 (1) on the real 
axis. When a Floquet multiplier A decreases through —1, 



the periodic orbit loses its stability via period-doubling 
bifurcation. On the other hand, when a Floquet mul- 
tiplier A increases through 1, it becomes unstable via 
pitchfork or saddle-node bifurcation. For each case of the 
period-doubling and pitchfork bifurcations, two types of 
supercritical and subcritical bifurcations occur. For more 
details on bifurcations, refer to Ref. [gj. 

Lyapunov exponents of an orbit {zm} in the Poincare 
map P characterize the mean exponential rate of diver- 
gence of nearby orbits |^^. There exist two Lyapunov 
exponents ui and (T2 {ai > a-i) such that cti -I- (T2 = — F, 
because the linearized Poincare map DP has a constant 
Jacobian determinant, detDP = e~^ . We choose an ini- 
tial perturbation bz^ to the initial orbit point zq and 
iterate the linearized map DP for Sz along the orbit to 
obtain the magnitute dm (= |^Zm|) of bz^n- Then, for 
almost all infinitesimally-small initial perturbations, we 
have the largest Lyapunov exponent oi given by 



(Ti = lim — In — — 

rn^oo m do 



(8) 



If (Ti is positive, then the orbit is called a chaotic orbit; 
otherwise, it is called a regular orbit. 



III. RICH DYNAMICAL BEHAVIORS OF THE 
DAMPED MO 

In this section, by varying the amplitude A, we inves- 
tigate the evolutions of both the stationary points and 
the rotational orbits of period 1 in the damped MO for a 
moderately damped case of F = 1.38. As A is increased, 
the damped MO, albeit simple looking, exhibits rich 
dynamical behaviors, such as symmetry-breaking pitch- 
fork bifurcations, period-doubling transitions to chaos, 
symmetry-restoring attractor-merging crises, and saddle- 
node bifurcations giving rise to new periodic attractors. 
In addition to these familiar behaviors, the stationary 
points also undergo a cascade of resurrections (i.e. an in- 
finite sequence of alternating restabilizations and desta- 
bilizations). It is found that the restabilizations (desta- 
bilizations) occur via alternating subcritical (supercriti- 
cal) period-doubling and pitchfork bifurcations. An infi- 
nite sequence of period-doubling bifurcations, leading to 
chaos, also follows each destabilization of the stationary 
points. 



A. Evolution of the stationary points 

We first consider the case of the stationary points. 
The damped MO has two stationary points f's. One 
is zi = (0,0), and the other one is zu = (|,0). These 
stationary points are symmetric ones with respect to the 
inversion symmetry 6*2, while they are asymmetric and 
conjugate ones with respect to the shift symmetry Si. 
Hence they are partially symmetric orbits with only the 



inversion symmetry 82- We also note that the two sta- 
tionary points are the fixed points of the Poincare map 
P [i.e., P{z) ^ z(z^ zj, zji)]. 

With increasing A we investigate the evolution of the 
two fixed points zj and £//. Two bifurcation diagrams 
starting from zj and zjj are given in Figs. |l|(a) and 
|l|(b), respectively. Each fixed point loses its stability via 
symmetry-conserving period-doubling bifurcation, giving 
rise to a stable 5'2-symmetric orbit with period 2. How- 
ever, as A is further increased each 5'2-symmetric orbit 
of period 2 becomes unstable by a symmetry-breaking 
pitchfork bifurcation, leading to the birth of a conju- 
gate pair of S'2-asymmetric orbits of period 2. (For 
the sake of convenience, only one S'2-asymmetric or- 
bit of period 2 is shown.) After breakdown of the 52 
symmetry, each 2-periodic orbit with completely bro- 
ken symmetries exhibits an infinite sequence of period- 
doubling bifurcations, ending at a finite critical point A* i 
(= 3.934 787 • • ■). The critical scaling behaviors near the 
critical point A* ^ are the same as those for the ID maps 
P, as will be seen in Sec. IV. 

After the period-doubling transition to chaos, four 
small chaotic attractors with completely broken symme- 
tries appear; they are related with respect to the two 
symmetries 5i and 52. As A is increased the differ- 
ent parts of each chaotic attractor coalease and form 
larger pieces. For example, two chaotic attractors with 
CTi (largest Lyapunov exponent) ~ 0.11, denoted by ci 
and C2, near the unstable stationary point zj are shown 
in Fig. H(a) for A — 3.937; their conjugate chaotic attrac- 
tors with respect to the 5i symmetry near the unstable 
stationary point i// are not shown. Each one is com- 
posed of four distinct pieces. However, as A is further 
increased these pieces also merge into two larger pieces. 
An example with cti ~ 0.18 is shown in Fig. 0(b) for 
A = 3.94. 

As A exceeds a critical value (~ 3.9484), the two 
chaotic attractors ci and C2 in Fig. 0(b) merge into a 
larger one c through an 52-symmetry-restoring attractor- 
merging crisis. For example, a chaotic attractor c with 
CTi ~ 0.37 and its conjugate one, denoted by s, with 
respect to the 5i symmetry are shown in Fig. |^(a) for 
A = 3.96. These two chaotic attractors c and s are 52- 
symmetric ones, although they are still 5i-asymmetric 
and conjugate ones. Thus the inversion symmetry 52 is 
first restored. However, as A increases through a second 
critical value (~ 3.9672), the two small chaotic attractors 
c and s also merge into a larger one via 5i -symmetry- 
restoring attractor-merging crisis, as shown in Fig. @(b) 
for A — 3.98. Note that the single large chaotic attractor 
with a 2± 0.64 is both Si- and 52-symmetric one. Con- 
sequently, the two symmetries 5i and 52 are completely 
restored, one by one through two successive symmetry- 
restoring attractor-merging crises. However, this large 
chaotic attractor disappears for A ~ 4.513, and then the 
system is asymptotically attracted to a stable rotational 
orbit of period 1 born through a saddle-node bifurcation, 
as shown in Fig. 0. 



B. Evolution of the rotational orbits 



We now investigate the evolution of the rotational or- 
bits of period 1. A pair of stable and unstable rotational 
orbits with period 1 is born for A ~ 2.771 via saddle-node 
bifurcation. In contrast to the stationary points, these 
rotational orbits are 5i-symmetric, but 52-asymmetric 
and conjugate, ones. The bifurcation diagram starting 
from a stable rotational orbit with positive angular ve- 
locity is shown in Fig. 0. (For convenience, the bifur- 
cation diagram starting from its 52-conjugate rotational 
orbit with negative angular velocity is omitted.) The 5i- 
symmetric rotational orbit of period 1 becomes unstable 
by a symmetry-breaking pitchfork bifurcation, which re- 
sults in the birth of a pair of 5i-asymmetric rotational 
orbits with period 1. (For the sake of conveninece, only 
one 5i-asymmetric orbit of period 1 is shown.) Then 
each rotational orbit with completely broken symme- 
tries undergoes an infinite sequence of period-doubling 
bifurcations, accumulating at a finite critical point A* 
{— 12.252 903 • • •). The critical behaviors near the accu- 
mulation point A* are also the same as those for the ID 
maps, as in the case of the stationary points. 

For A > A*, four chaotic attractors with completely 
broken symmetries appear; they are related with respect 
to the two symmetries Si and 52. Through a band- 
merging process, each chaotic attractor eventually be- 
comes composed of a single piece, as shown in Fig. ^(a) 
for A = 12.32. Four chaotic attractors with ai ~ 0.36 
are denoted by ci, C2, si, and S2, respectively. How- 
ever, as A passes through a critical value (~ 12.3424) the 
four small chaotic attractors merge into a larger one via 
symmetry-restoring attractor-merging crisis. An exam- 
ple for A = 12.38 is given in Fig. |(b). Note that the 
single large chaotic attractor with ai ~ 0.64 has both 
the 5i and 52 symmetries. Thus the two symmetries 
are completely restored through one symmetry-restoring 
attractor-merging crisis, which is in contrast to the case 
of the stationary points. 

However, as shown in Fig. 0, the large symmetric 
chaotic attractor in Fig. p|(b) also disappears for A ~ 
13.723, at which saddle-node bifurcations occur. After 
disappearence of this large chaotic attractor, the sys- 
tem is asymptotically attracted to a stable 52-symmetric, 
but 5i-asymmetric, orbit of period 2 born via saddle- 
node bifurcation. This stable 52-symmetric orbit with 
period 2 also exhibits rich dynamical behaviors simi- 
lar to those of the stationary points. That is, as A 
is increased, a symmetry-breaking pitchfork bifurcation, 
period-doubling transition to chaos, merging of small 
asymmetric chaotic attractors into a large symmetric one 
via symmetry-restoring attractor-merging crisis, and so 
on are found. However, unlike the cases of the station- 
ary points and the rotational orbits, the large symmetric 
chaotic attractor disappears for A = 23.751799 •••, at 



which the two unstable stationary points zj and zu be- 
come restabihzed through subcritical period-doubhng bi- 
furcations. These "resurrections" of the stationary points 
wiU be described below in some details. 



C. Resurrections of the stationary points 

The linear stability of the two stationary points zj and 
zji is determined by their linearized equations, 



Sx + r SxT 27ryl cos 2TTt Sx — 0, 



(9) 



where — {+) sign of the third term corresponds to the 
case of zj (zn). (The linearized equation of zjj can also 
be transformed into that of zj by just making a shift in 
time, t — > i -I- i.) Note that Eq. (Q) is just a simple form 
of the more general damped Mathieu equation [|lj] . It is 
well known that the Mathieu equation has an infinity of 
alternating stable and unstable A ranges. Hence, as A is 
increased, the stationary points undergoes a cascade of 
"resurrections," i.e., they will restabilize after they lose 
their stability, destabilize again, and so forth ad infini- 
tum. 

It is found that their restabilizations (destabiliza- 
tions) occur through alternating subcritical (supercriti- 
cal) period-doubling and pitchfork bifurcations. As ex- 
amples, we consider the first and second resurrections 
of the stationary points. The first resurrection of i/ is 
shown in Fig. ||. When A passes through the first resta- 
bilization value {— 23.751799 • • •), the rightmost large 
symmetric chaotic attractor in Fig. M disappears and the 
unstable stationary point zj restabilizes via subcritical 
period-doubling bifurcation, giving rise to an unstable or- 
bit of period 2. Two bifurcation diagrams starting from 
the restabilized zj and zu are given in Fig. H(a) andy(b), 
respectively. Each stationary point loses its stability via 
symmetry-breaking pitchfork bifurcation, giving rise to 
a pair of S'2-asymmetric orbits with period 1; only one 
asymmetric 1-periodic orbit is sho wn. This is in con- 
trast to the case given in Sec. [II A (compare Fig. |]with 
Fig. P) , where the stationary points become unstable via 
symmetry-conserving period-doubling bifurcations. Af- 
ter breakdown of the ^2 symmetry, an infinite sequence 
of period-doubling bifurcations follows and ends at its 
accumulation point A* 2 (= 24.148 001 • • •). When A ex- 
ceeds A* 2 , a second period-doubling transition to chaos 
occurs. The critical scaling behaviors of period doublings 
near the second critical point A = A^ 2 are also the same 
as those near the first critical point A* i . 

Dynamical behaviors of the damped MO after the 
second per iod-doubling transition to chaos are shown 
in Fig. |lO|(a). As A passes through a critical value 
(~ 24.1549), small chaotic attractors with completely 
broken symmetries merge into a large symmetric chaotic 
attractor via symmetry-restoring attractor-merging cri- 
sis. However, the large chaotic attractor also disappears 



for A ~ 29.342, at which saddle-node bifurcations oc- 
cur. After disappearence of the large chaotic attractor, 
the damped MO is asymptotically attracted to a stable 
S'l-symmetric, but 5'2-asymmetric, orbit of period 1 born 
via saddle-node bifrucation. The subsequent evolution of 
the stable 1-periodic orbit is shown in Fig. ^(b). Note 
that it is similar to that of the rotational orbit described 
in IIIB (compare Fig. [lO|(b) with Fig. 0). 

The rightmost large symmetric chaotic attractor in 
Fig. |l^(b) appears via symmetry-restoring attractor- 
merging crisis for A = 57.67. However, it also disap- 
pears for A = 67.076 913 • • •, at which each stationary 
point restabilizes again. Unlike the case of the first res- 
urrection (see Fig. ||), this second resurrection of each 
stationary point occurs via subcritical pitchfork bifurca- 
tion, giving rise to a pair of unstable 1-periodic orbits 
with broken symmetries. The second resurrections of the 
two stationary points zj and z// and their subsequent 
bifurcation diagrams are shown in Fig. |ll|(a) and |ll|(b), 
respectively. Note that these third bifurcation diagrams 
are similar to those in Fig. ||. The critical scaling behav- 
iors near the third period-doubling transition point A* 3 
{— 67.104 872 • • ■) are also the same as those near the 
first period-doubling transition point A = A* ■^. 



IV. CRITICAL SCALING BEHAVIORS IN THE 
PERIOD-DOUBLING CASCADES 

In this section, we study the critical scaling behaviors 
in the period-doubling cascades. The orbital scaling be- 
havior and the power spectra of the periodic orbits born 
via period-doubling bifurcations as well as the parameter 
scaling behavior are particularly investigated. 

The critical scaling behaviors for all cases studied are 
found to be the same as those for the ID maps. As 
an example, we consider the first period-doubling tran- 
sition to chaos for the case of the stationary points. As 
explained in Sec. Ill A, each stationary point becomes 



unstable through symmetry-conserving period-doubling 
bifurcation, giving rise to a stable S'2-symmetric orbit of 
period 2. However, each S'2-symmetric orbit of period 
2 also becomes unstable via symmetry-breaking pitch- 
fork bifurcation, which results in the birth of a conju- 
gate pair of S'2-asymmetric orbits with period 2. Then, 
each 2-periodic orbit with completely broken symmetries 
undergoes an infinite sequence of period-doubling bifur- 
cations, ending at its accumulation point A^ i. Table H 
gives the A values at which the period-doubling bifurca- 
tions occur; at Ak, a Floquet multiplier of an asymmetric 
orbit with period 2'' becomes —1. The sequence of A^ 
converges geometrically to its limit value A* ^ with an 
asymptotic ratio S: 

4 ^ 4^_li^ ^ S, (10) 



A 



fe+i 



A, 



The sequence of 6k is also listed in Table I. Note that its 
limit value S (~ 4.67) agrees well with that (= 4.669 • • •) 



for the ID maps . We also obtain the value of A* ^ 
(= 3.934 787024) by superconverging the sequence of 

As in the ID maps, we are interested in the orbital 
scaling behavior near the most rarified region. Hence, we 
first locate the most rarified region by choosing an orbit 
point z^'^^ [= (j:^'^-', y'-'^^)] which has the largest distance 
from its nearest orbit point P^ {z^'^'') for A = A^. The 
two sequences {x^^^} and {y*-*^-*}, listed in Table O, con- 
verge geometrically to their limit values x* and y* with 
the ID asymptotic ratio a (= —2.502 • ■ •), respectively: 



Jk) 



n(fc-l) 



,(fc) 



,,(fc-l) 



dx.k 



Jk+i) 



Jk) 



a, Oiy,k 



,(fc+i) 



.Ak) 



(11) 

The values of x* (= 0.015 369) and y* (= 0.348 175) are 
also obtained by superconverging the sequences of x'^'^^ 
and y^^'^ , respectively. 

We also study the power spectra of the 2''-periodic or- 
bits at the period-doubling bifurcation points Ak- Con- 



sider the orbit of level k whose period is g = 2*^, {z 

{xm\ym), m = 0, 1, ... ,(7 — 1}. Then its Fourier 
ponent of this 2'^-periodic orbit is given by 



.(fc) 



K) = - 



1^ 



Jk) -iujjm 



(12) 



m=0 



where ojj = 2'kj /q, and j = 0, 1, . . . , g — 1. The power 
spectrum P^^^{ujj) of level k defined by 

P«(c.,) = |zW(c.,)p, (13) 

has discrete peaks at lu — ujj. In the power spectrum 
of the next (fc + 1) level, new peaks of the (k + l)th 
generation appear at odd harmonics of the fundamental 
frequency, luj = 27r(2j + l)/2('=+i) (j = 0, . . . , 2'^' - 1). To 
classify the contributions of successive period-doubling 
bifurcations in the power spectrum of level A;, we write 



P('^-) - Poo5{lo) + Y, 



k 2('-i)- 

E 

1=1 j=a 






p>:'5{uo 



UJlj] 



(14) 



^(fe) 



where PA- is the height of the jth peak of the /th gen- 
eration appearing at w = ujij (= 27r(2j -I- l)/2'). As an 
example, we consider the power spectrum P^^> [lo) of level 
8 shown in Fig. ^ The average height of the peaks of 
the /th generation is given by 



</>W(0 = 



1 



2('^ 



2'--'-! 

_ V^ p(fe) 

1) Z^ ^i] ■ 

J=0 



(15) 



It is of interest whether or not the sequence of the ratios 
of the successive average heights 



?('=)(';^ - ^C^) 



2(3^''>{l)^q^^''>{l)/<j)^''>{l + l), 



{k)l 



(16) 



converges. The ratios are listed in Table [II. They seem 
to approach a limit value, 2/3 ~ 21, which also agrees well 
with that (= 20.96 • • •) for the ID maps |1. 



V. SUMMARY 



Dynamical behaviors of the damped MO are investi- 
gated in details by varying the amplitude A. The damped 
MO, despite its apparent simplicity, exhibits rich dynam- 
ical behaviors such as breakdown of symmetries via pitch- 
fork bifurcations, period-doubling transitions to chaos, 
restoration of symmetries via attractor-merging crises, 
and creation of new periodic attractors through saddle- 
node bifurcations. In addition to these familiar behav- 
iors, a cascade of resurrections (i.e., an infinite sequence 
of alternating restabilizations and destabilizations) of the 
stationary points occurs. It is found that the stationary 
points restabilize (destabilize) via alternating subcriti- 
cal (supercritical) period-doubling and pitchfork bifurca- 
tions. An infinite sequence of period-doubling bifurca- 
tions, leading to chaos, also follows each destabilization 
of the stationary points. The critical scaling behaviors 
in the period-doubling cascades are found to be the same 
as those of the ID maps. 
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TABLE L Asymptotically geometric convergence of the 
parameter sequence {Ak} 



Sk 



1 3.911404100 371 

2 3.929 795 227 873 4.690 

3 3.933 716 964 019 4.664 

4 3.934 557 747 089 4.667 

5 3.934 737 918 915 4.669 

6 3.934 776 506 700 4.668 

7 3.934 784 773 689 4.673 

8 3.934 786 542 923 



TABLE II. Asymptotically geometric convergence of the 
orbital sequences {x^^'} and {y }■ 



k 


x^"^ 


ax,k 


,(.) 


ay,k 


1 


0.012 394 993 




0.337 125 704 




2 


0.016 703 927 


-2.354 


0.350 786 135 


-3.410 


3 


0.014 873 276 


-2.607 


0.346 780 101 


-2.133 


4 


0.015 575 414 


-2.451 


0.348 658 607 


-2.714 


5 


0.015 288 893 


-2.532 


0.347 966 417 


-2.395 


6 


0.015 402 057 


-2.487 


0.348 255 392 


-2.561 


7 


0.015 356 562 


-2.511 


0.348 142 575 


-2.472 


8 


0.015 374 678 




0.348188 212 





TABLE m. Sequence 2l3'-'''>{l) [= (t>'-^\l) / (l>'-^\l + 1)] of 
the ratios of the successive average heights of the peaks in 
the power spectra 
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FIG. 1. Bifurcation diagrams starting from (a) the 
S2-symmetric, but Si-asymmetric, stationary point zj and 
(b) its Si-conjugate stationary point zn. The first and sec- 
ond P2's denote the stable ^-ranges of the 5'2-symmetric and 
S'2-asymmetric orbits of period 2, respectively. The other 
PN (A'^ = 4, 8) also designates the stable A-range of the 
Sg-asymmetric periodic orbit with period N. 



FIG. 2. Band-merging of chaotic attractors. (a) For 
A — 3.937, each of the chaotic attractors Ci and C2 is com- 
posed of four pieces. However, as A is increased these pieces 
also merge to form two larger pieces. An example for A = 3.94 
is shown in (b). 
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FIG. 3. Attractor-merging crises. The two chaotic attrac- 
tors ci and C2 in Fig. 2(b) merge into a larger one c via 
S'2-syninietry-restoring crisis. For A = 3.96 the chaotic at- 
tractor c with the inversion symmetry S2 and its conjugate 
one s with respect to the Si symmetry are shown in (a) . These 
two Si-asymmetric small chaotic attractors c and s also merge 
to form a larger one via Si -symmetry-restoring crisis. A sin- 
gle large chaotic attractor with completely restored Si and 
S2 symmetries is shown in (b) for A = 3.98. 
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FIG. 4. Jump to a rotational orbit. The large symmet- 
ric chaotic attractor in Fig. 3(b) disappears for A ~ 4.513, 
and then the asymptotic state of the damped MO becomes 
a stable rotational orbit with period 1 born via saddle-node 
bifurcation. Such a saddle-node bifurcation, giving rise to 
a pair of stable and unstable orbits of period 1, occurs for 
A ~ 2.771 (a stable one is denoted by a solid line, while an 
unstable one is represented by a dashed line). A bifurcation 
diagram starting from zi is also shown. 
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FIG. 5. Bifurcation diagram starting 

from the Si-symmetric, but 5'2-asymmetric, rotational orbit 
with period 1. The first and second Pi's denote the stable 
^-ranges of the Si-symmetric and 5'i-asymmetric orbits of 
period 1, respectively. The other PN (A'^ = 2,4, 8) also desig- 
nates the stable A-range of the Si -asymmetric periodic orbit 
with period N. 
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FIG. 6. Attractor-merging crisis for the rotational case. 
Through a band-merging process, each of the four chaotic 
attractors c\, C2, si, and S2 with completely broken symme- 
tries eventually becomes composed of a single piece, as shown 
in (a) for A — 12.32. The four small asymmetric chaotic at- 
tractors merge into a larger one via symmetry-restoring crisis. 
A single large chaotic attractor with simultaneously restored 
Si and S2 symmetries is shown in (b) for A = 12.38. 
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FIG. 7. Jump to an orbit with period 2. The large symmet- 
ric chaotic attractor in Fig. 6(b) disappears for A ~ 13.723, at 
which a saddle-node bifurcation, giving rise to a pair of stable 
and unstable orbits with period 2, occurs. (A stable one is 
denoted by a solid line, while an unstable one is represented 
by a dashed line) . Then the damped MO is asymptotically at- 
tracted to the stable S'2-symmetric, but Si-asymmetric, orbit 
with period 2. This stable 5'2-syininetric orbit also exhibits 
rich dynamical behaviors similar to those of the stationary 
points. In the left part, a bifurcation diagram starting from 
a rotational orbit with period 1 is also given. 



FIG. 8. First resurrection of the stationary point zj. When 
A passes through the first restabilization value (~ 23.752), the 
unstable stationary point zj restabilizes through a subcritical 
period-doubling bifurcation, giving rise to an unstable orbit 
with period q = 2. The solid and dashed lines also denote 
stable and unstable orbits, respectively. 
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FIG. 9. Seoncd bifurcation diagrams starting from (a) 
the restabilized stationary point zi and (b) its 5'i-conjugate 
stationary point in. Here the PN designates the stable 
^-range of the 5'2-asymmetric periodic orbit with period N 

(N=l,2,4,8). 



FIG. 10. Dynamical behaviors after the second pe- 
riod-doubling transition to chaos is shown in (a). For 
A ~ 24.1549 small chaotic attractors with completely broken 
symmetries merge into a large symmetric chaotic attractor 
via symmetry-restoring crisis. However, this large symmet- 
ric chaotic attractor disappears for A ~ 29.342, and then 
the damped MO is asymptotically attracted to a stable or- 
bit of period 1 born via saddle-node bifurcation. As shown 
in (b), subsequent evolution of the stable Si-symmetric, but 
S'2-asymmetric, 1-periodic orbit is similar to that of the rota- 
tional orbit shown in Fig. 7. Here the solid and dashed lines 
also denote stable and unstable orbits, respectively. For other 
details see text. 
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FIG. 12. Power spectrum P'^^\u}) of level 8 for ^ = As 
(=3.934 786 542 923). 



FIG. 11. Second resurrections of the stationary points and 
third bifurcation diagrams starting from (a) the restabilized 
stationary point zi and (b) its Si-conjugate stationary point 
zii. When A passes through the second restabiliztion value 
(~ 67.08), each unstable stationary point becomes restabi- 
lized via subcritical pitchfork bifurcation, giving rise to a pair 
of unstable orbits with period g = 1. The solid and dashed 
lines represent stable and unstable orbits, respectively. The 
third bifurcation diagrams are similar to those in Fig. 1; the 
symbols are also the same as those of Fig. 1. 
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